home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / DSIFA.z / DSIFA
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDSSSSIIIIFFFFAAAA((((3333FFFF))))                                                            DDDDSSSSIIIIFFFFAAAA((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DSIFA   - DSIFA factors a double precision symmetric matrix by
  10.      elimination with symmetric pivoting.
  11.  
  12.      To solve  A*X = B , follow DSIFA by DSISL.  To compute  INVERSE(A)*C ,
  13.      follow DSIFA by DSISL.  To compute  DETERMINANT(A) , follow DSIFA by
  14.      DSIDI.  To compute  INERTIA(A) , follow DSIFA by DSIDI.  To compute
  15.      INVERSE(A) , follow DSIFA by DSIDI.
  16.  
  17.  
  18. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  19.       SUBROUTINE DSIFA(A,LDA,N,KPVT,INFO)
  20.  
  21. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  22.      On Entry
  23.  
  24.      AAAA DOUBLE PRECISION(LDA,N)
  25.         the symmetric matrix to be factored.
  26.         Only the diagonal and upper triangle are used.
  27.  
  28.      LLLLDDDDAAAA INTEGER
  29.         the leading dimension of the array  A .
  30.  
  31.      NNNN INTEGER
  32.         the order of the matrix  A .  On Return
  33.  
  34.      AAAA a block diagonal matrix and the multipliers which
  35.         were used to obtain it.
  36.         The factorization can be written  A = U*D*TRANS(U)
  37.         where  U  is a product of permutation and unit
  38.         upper triangular matrices, TRANS(U) is the
  39.         transpose of  U , and  D  is block diagonal
  40.         with 1 by 1 and 2 by 2 blocks.
  41.  
  42.      KKKKPPPPVVVVTTTT INTEGER(N)
  43.         an integer vector of pivot indices.
  44.  
  45.      IIIINNNNFFFFOOOO INTEGER
  46.         = 0  normal value.
  47.         = K  if the K-th pivot block is singular.  This is
  48.         not an error condition for this subroutine,
  49.         but it does indicate that DSISL or DSIDI may
  50.         divide by zero if called.  LINPACK.  This version dated 08/14/78 .
  51.      James Bunch, Univ. Calif. San Diego, Argonne Nat. Lab.  Subroutines and
  52.      Functions BLAS DAXPY,DSWAP,IDAMAX Fortran DABS,DMAX1,DSQRT
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.